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THE CALCHLAIDION 01? THE THKItMAL UTILIZATION 
IN CELIfi WITH COMPLEX FUEL ELEMENTS 


G*A,Bat, E, A, Grigoryeva, V,N.Iiobedev , A»l. Prokhorov 
L.S, Tzygajikov 

I. 

It seems possible that major errors in the calculation 
oi‘ 6- coefficient of the thermal utilisation in the hetero- 
geneous nuclear reactor lattices should be determinedly- 
excluded by means of a transition from a ane-velocity diffu- 
sion equation to a klneti'- "nuation of thermal izat ion with 
its rather detailed description of the cross-sections of the 
introduction between neutrons and chemically bound moderator 
nuclei in a thermal movement. The research dedicated to this 
programm has been going on for already a few years ant, 
in fact, are nearing completion now. 

Yet, in the cases of numerous calculation versions, 
those of nuclear fuel burnout included, the overall considera- 
tion of kinetic and thermalizatlon effects is extremely 
oumberaome. This primarily concerns lattices containing not 
only separate cylindrical fuel elements but rather whole 
bundles of them or a system of co-axial tubes, in other words, 
channels, surrounded by a moderator layer which is thicker 
than thermal neutron free path. This chaumel structure is 
typical to the Lenin icebreaker reactor and the reactor MR, 

The consideration of the channel's geometry and changing of 

its composition in the course of the burn-out needs calculati?< 
programmes whose realisation takes a good many hours of fast 

computer operation, She necessity then springs up to evolve 
sipmler methods to ensure satlsfacory precision and to rest 
on reference data obtained in experiments and accurate cal- 
culating procedures. It looks feasible for this purpose to 
tackle -tire problem in two stages: to determine fii*st of all 
the space variations of the thermal neutron spectrum form or 
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eleo, which In Ifie equivalent, the space dependence of cross- 
eectloiis then take up the one-velocity equation with the 
pre-civea Law of average croBs-sectionc changeOi furthermore, 
it makoe lota of eenao to try the siapieet model to calcula- 
te the neutron spectrum, for example, the diffusion heavy- 
-gas model (porhapa, with effective moderator parameters), 
and seek to find an acceptable approximate solution of the 
thermalization equation. Bearing it in mind that the thormall- 
zatlon correction is commonly as to the value (1-0 )// 

we decreaee an error by an order in the case of sys- 

tems without high temperature gradients, free of substances 
with low energy resonances. The problem's second peirt for comp- 
licated lattices, one cannot obviously do without the kinetic 
equation to be solved, at least one per a cell, and two 
or three times at a deep burnout. Thus , the artificial 
means of the typo of the Amouyal — Benoist-Gorowltz integral 
method &r 9 practically out of the question. Expansion into 
spherical harmonics does not necessarily provide the duo accu- 
racy, the Monte— Carlo methods are rational only to low symmet- 
ry charnels (for those, say, of the first atomic power sta- 
tion). The characteristics method (Vladimirov) and Sn-methods 
(Oarlason) seem most suitable for symmetric cells. 

2 . 

The equation of neutron thermalization in the reactor 
coll moderator in the space of coordinates and velocities un- 
der diffusion approximation (it is supposed the mmoatomic 
heavy gas model) assumes the form: 


DAf ('",'') + ■ 



[^Is 


2 


kTcp .,3 ()f , ^,^3 .,4 


m '' 


+ 


vM 


rvT;f^Y^Qo(v)=u. 


( 1 ) 

Here -f (^,Y) is the neutron density in the space of coor- 
dinates and velocities ) 

D=(3Z^p)*'^ ordinary macroscopic con- 

stants, which a 3 ?o said to bo independent of t!ie velocity V \ 
- temperature of a moderator in °K, K - the constant 

of Soltzaan. a - sass of neutron, ^ -average logary thmic 
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energy decre'uenU per coliision; - Uio moot probable velo- 

city for the Majcweliian opoctruixj o" bho temperature T ; 

- noutroa soiirco. 

Lot U8 muitipiy equation (I) by ■/' i!v oud integrate for 
all the volocitiee . Aware of thr' h';iiaA'ior of -f at M =0 
and oo , obtain: 

N(r) + (^--0, (2) 

oo ^ 

where N (?) = | V^t (r.V ) d V, 9(r) - | V^f (r, V) dv - are the 

O o 

integral density and the flujc of neutrons. 

Equation (2) differs from an ordinary one- velocity dif- 
fusion equation for thermal neutrons in its second member | 
which substitutes for the member ( vVv )Z!^ (r ) . There- 

fore, a necessary suggestion for putting down a one-velocity 
diffusion equation reads that the neutron integral density 
h/fr) can be substituted by the ratio ^(r)/V , where 

V - the average velocity for the Majcwellian spectrum of a 
certain temperature Tpg independent of the coordinates 
and ascribed to the neutron distribution. It will be only na- 
tural to make here a next step on the suggestion still there 
about, the Maxwellian type of the neutron spectrum in the 
cell's moderator^ the spectrum temperature being looked upon 
as a coordinate function. It should be noted that the results 
of calculations and eacperimental data both confirm this appr 
roxlmation. 

By the multiplication of equation (I) by V^dv and 
integration on all the velocities it is obtained: 

where H^(P) = | f (r, v) dv. 

o 

Substituting in equation (2) , (3) ‘^'(?)/v(r) for 

liCr) amd taking it Into account tiiat the spectrum nature Cica of 
neutrons is similar in all the moderator points, the system 
emerges 
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DAr(r)-^<p + ^.^- (5) 


-8 l±.f-± 'tn.s. 

'3ar 9(r) ’ I v°i tt 

The relation between the free members q and is not 
difficult to get making use, for one, of the weil-known result 
pertaining to the spectra in homogeneous mediums without apa- 
ce gradients s 

K=i,84 

and the corresponding solution of the system (4) + (6) 

’^n.s. - * I (^ " ' ^] 

This brings = 2,^08 qV®, 

The boundsury conditions on the cells* Sp outer surface 
are apparent! 

Once within the fuel block there are no light nuclei and the 
exchange of energT' between neutrons and nuclei is non— existent , 
the purely diffusion boundau^ conditions on the surface of the 
cylindrical block section and moderator can read as follows: 

^V<Pl J tn f-fDo IJ 


Where Do=(2>2or , U “ = ^oc (v) + los -usually ma- 

croscopic fuel parameters. 

Conditions (9) and (10) are easily generalised as to other 
fuel block geometries. 

The above calculation pro<5edure for a neutron thermaliza- 
tion model on the ideal monofttomic gas of heavy nuclei can be 
applied to real mediums carrying the selection of the effecti- 
ve values of K and H parameters* 

Taking into account a very peculiar calculation of the 
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water moderator it must be bomo ! - mind not only how rough 
the heavy gas model la, but alno i;r:« auggestion madei D = 
coasts The correspondlxig account o. dependence D(v) would 

cause the change of the equation ('O, ( 5 ) without giving them 
the due accuracy. Proceeding from -he asaumption that in the 
water medium D -^v and introducing immediately into the app- 
roximate equations two constants ejid = D 

to be considered independent of the coordinates, it is found: 

Dh. [Tee«’dK]{rx‘^e-‘“clx] 

D<p [TxV'dxJ' 

The choice of effective K and D values will be given in 
section 3 . The aproximate solution obtained above of a therma- 
lization equation makes it possible to average the cross-sec- 
tions of neutron interaction with the cell substances over the 
and 

energy range then to solve the one-velocity kinetic problem 
with the cross-sections dependent on coordinates. It is note- 
worthy, that there is no implication of the Maxwellian spect- 
rum form witMn the block whatsoever. 

The one-velocity kinetic equation for the neutron distri- 
bution function in the cell p is given the form: 

^'7F’-!:F-^fF(r,8')9(JJo)d9'+S(r). ( 11 ) 

Here r is the vector-radius of points* x, 'j, z; Q 
- a single vector of the neutron velocity direction,^Q=cos(Si-9') ; 

g(jij is the indioatxlx of scattering , S(r) “is 
the source of neutrons. 

As usual, Q=i 2 ,(+j 9 y+kS 22 *, 9x=sin\3’'COS9i Qy=sin\5'sin(p; S? 2 =cos I?", 

where -is the angle between K and Q , while -is 
the angle between I and the projection Q against the 
plane x, ^ 

To solve the equation (II) the method of characteristics 
has been used fs) and the code has been compiled 

for a fast computer. 

Straight line segments - the paths of neutron constitute 
the characteristics of the linear differential expression QVF . 
Then the straight line I is chosen in coordinate space, and 
* 1 ^ - parameter is introduced descriptive of the length on 1 . 

'he points of the straight line the neutron flux moving in 
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tjie direct;! on , coinii.t JLing with 

in terms of the equationi 

^■HF=^[F(r,C')9(Po)d»VS(f) . 

For obtaining the exhausting picture of neutron distribution 
in the circular cylindrical cell it is sufficient to consider 
the straight lines 1 on the planes y = const, and neutrons 
moving along those straight lines. Indeed, the function P in 
this case is dependent on the module of projection r against 
the plane x, y, , the angle 0 between the direction of neutron 
movement and the axis z eind the angle between projections l 
and 5? against the plane x, y. Apparently, the variation range 
of these variables is exhausted as long as we move along all 
the straight lines 1, 

■tet, the characteristics 1 forms the angle © with 
the axis z, then = (1 - where ^ = coaO , and 

the equation (12) may be written as follows 


I 


is describabl* 


"57 






( 13 ) 


The neutron reflection conditions on the cell boundary 
are looked upon as the boundary conditions for P 

is the cell radius, 

ja=x/r = cos<P; 0« 

It is possible to consider that cross-sections 
step constant functions, 

and the scattering is linear anisotropic one that is 

g(Mo) = l + 3jio Mo, 

where jIq represents the average scattering angle cosine. 
The right part of the equation being expressed (|3) through Q, 
the equation reads as 
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Is9 + 3}ioM CI+ S(r') 
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The equation ('. !•; with the 
solve by Lho Iteration meUi(A.. 
xiiuate expretJBion for ■ f 


tion contaixTis P, 


1 


value o.r 


'll) ia convenient to 
e prodotermined appro- 
t,do first approxiiaa- 

to come 
-so- 

:e.li cnIcnJation. The 


in and ’J , 
. ut(.t I'll tu re to 


over, however, in the moot v v 

lubiori aB the fluai one ..or 

ordinary nonhowr^onius differential equabioa of the first 
order along the characteriatic I 

^ L F ( K 'I Y 'l = 0 

is solved at the fixed ^ and y parameter values, P being 
determined in the cross-points between the choi’acteristic and 
t;he straightline X=0 and somi-circumf erences ^ 


The variable ^ 


mesh-points find themselves to be 
the positive roots of the Legandre even polynomials, mean- 
while, the spacing on j is not affected by the choice • 

The solution of the equation (16) calls, above all., for 
finding P=Pq a-t the characteristic's left end. For this pur- 
pose, in the equalities, which link up the two function va- 
lues F in meighbouring points on the characteris.tlc, all 
the are subsequently excluded from as early as F^. As a re- 


and comes about where 

n n 


suit, the relation between P^ 
the value of F at the right end. The relation, as well as 
the boundary condition Fq=Fj| defines P^. 

For the distribution function integral the formula is 
taken of the fourth precision order* 

Once integrated, the equation (16) along the characte- 
ristic within the stretch from point ( , )j ) to point 


is 


> y ) yields 


F(>:n.V4)- F(Xn.pV.^j) + ^7^ j = | Qdx , 


T 


where 


Xp ..| 


Fa.=^h4 



^n-1 Xn-l 

^)+0fh*) 


( 17 ) 
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On solving the equation (16) with bhe right part the cycle 
resumes as soon as the new value of Q is arrived at. To bring 
about better convergence the iteration method is used whore 
the major deviation from the exact solution is to be found 
within the - approximation, and the L.A. Lusternik method 
is involved as supplementary one* C3l. 

The process is ended ae soon us the difference of va- 
lues in ooy of the cell points for the two forthcoming itera- 
tions comes below the given value £ . ^ being equal to 
10 it is sufficient accomplish ~ 15 iterations. The 

functions F, and D come to the print. It takes all 

in all 10 - 15 minutes to calculate the cell. 

Up to 12B points along r , 80 point for ^ and 16 
points for j are allowed. 

In further calculations it proved indispensable to use 
as many as 100 points along the radius, 60 y -points and 5 
points for ^ • The vicinity of the moderator and absor- 

ber section boundaries, that Is to say, in the areas of sub- 
stantial changes of T , the radial step was smaller. 

The s^>aclng for y grows in the absorber, because of the 
h- ’fir flux anisotropy. 

It IS noteworthy, that the specific one-velocity 
calculation of a few ring multi-layer cells with a water 
moderator- the cell study is forth-coming, - was indicative 
of that the accepted neutron scatterl.g anisotropy count was 
detailed in excess, and is liable to an easy subs- 
titution by = (1" )lo)^S without any loss of accuracy 

as is done in the diffusion approximation. 

As the reference point in controlling the code the 
calculation of a double-zone cell was chosen. Investigated 
in the stiidy 

Its parameters:' Ro = 1i3cm, R. = 11,28 cm , 

Ioc=0,3230 cmi Eo=0,722lcm"^ = 

“0,3721 cm"*’, Jioo=jloi =*0 , So=0^ S^ = l 

The below result of the calculation of the characteristic 
method ratio is compared with the outcomes of the 

study (4) (Table I ) 

It fflUBh be borne in mind that; the reflection conditions 
on the cell boundary (14) are not rigorous, provided the 

water film ^ tMn enough the latter might involve higher 
values for W/% , 

V . A — . 
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In the above Inatanco this ercor does not ahow throu^, 
jet, tlia reaulta of the folio* in,/ oection for the aaall radluB 
cells are likely be pregoant ihe lack of perelsioa. 

The thexsal utlliaationi caJoulj»iliou scheae as described in 
section 2 requires, as aentioned above, the selection of ef- 
fective values for theiaalieation paraaeters K and D. 

Fig, I shove the calculation results of neutron gas 
teaperatux^ distribution in the cell*^ with a water Moderator 
and block of the natural uraniiua, detailed upon in the study 
bj Mostovoy The teaperature of the intra-biock neutron gas 
under the conditions of a stMWinaeawlliaa fora spectrvia was 
understood to be the value of '^'n.y where V[p 

the Most probable velocity of the distribution h (v) • 

It is coaaon knowledge that the experimentally obtained 
neutron spectrun was that of neutrons aoving parallel to the 
blocks; therefore as long as r = 0 the stars aark down the 
values obtained by Martchuck, Siaelov (etc.)^°J in solving the 
kinetic theraallzation equation for parallel and normal 
beaus and for spectrum which hajs been averaged in the direc- 
tion of the neutron flight. It is evident, that the values 
K = 1,46 (7) /respectively , q^ =2,223 q / 

and D = 0,1561 ^Tn,g,/29^ ensure a substantial degree 

of coincidence of our results wi th the mort reliable infor- 
aation pertaining to the neutron s^ectrua in this cell* 
Instaintaneously upon the finding of (for the 

indoor water teaperature only) followed the determination of 
the N(r) - integral density distribution represented in 
Fig. 2 by the dotted line. The aolid line there is a result 
of the kinetic calculation with the prematurely selected 
distribution ^ (r); the points are the results of the 

experiments by Stolyarov and al. /8/ obtained by way of cop- 
per foil activation. Conscious of the experimental errors 
outside the drawing, the calculation still looks satisfactory; 
the diffusion curve proved too sloping as should have been 
expected* 


a calculation was made of a cell with the set aJ-uminina 
Che exact distribution of alluminium being represented 

i73 
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The "exact” value of 6 = 0 , 7 ^ the diffusion calculation 
gives 0 = 0,BOO7 and, o^auecjuently , brin^o about an 

error 10^ into the ratio 


The Bchemo also underwent testing by the calculation of 
other geometrically simple uran lum-water lattices, among the* 
those with slightly enriched metallic blocks which once were 
a matter of most close study in the Brookhaven National la- 
boratory, USA (BNL)* Tables II /md III quote the typical re- 
sults for the cases where the "blackness” is nearing the 
extreme* 

It should be remembered that throughout our calculus the 
use was made of the uranium neutron cross-section system fio], 
slightly different from that accepted in the Brookhaven labo- 
ratory, all the more since the transition to the latter con- 
siderably deteriorates the agreement between the experimental 
and calculation datse It is very likely to be explained by 
the insesitivity of our calculations toward the position of 
the upper boundary of the thermal neutron spectrum* 

For simple cells "the pre-heating" within the block is 
practically independent of the lattice step, and the tempe- 
rature of a cell boundary is approximately desdribed by the 
0laple relation r , 

It is of interest to trace whuitev.? 'effect different calcu- 
lation approximations ]invo on the caiw od value of 0. Table IT 
shows an example of t-hlr^ on tin-. • ' of cella with 

, • clrp, of 0*600 inch diamijtor r.n.i 1.3% , (See 

Giearly, in all the above cases the diffusion approxiaa- 
tion enror ia although tbers^lization effect at 

t;ho large spaeixu^ over greater l??i>ortonce • 

Thus, cyllndric!^} hiock calculations confirm 

efficiency :f this x:,otl 70 d i r.v- th€7 door open for consi- 
t : * rr aiortj coapxicated calls* oncrpl.e, the light-water 

been chOi:^oxi with fuel channels asp?:aially designed 
.^viidy of '/ tiir reactor cjxi tlcallty* Tlio chaxt- 

j j r/hjole of the 


• rbn 


- iijilng sjptea of two, 
■ being ijiseirted 
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into each other. 

UOj with the enrichaent of ~ 30% nerve ae fuel, the clad- 
ding and the channel 'a jacket tube are aade out of otalnleaa 
steel, the aizes transpire very eoll froa the figures. 

fig. 3. and 4 illustrate the results of an one-veioclty 
calculation for the flujt distribution along the three— ring 
channel ceil, as effected in ^ 5 “ and kinetic appro- 

xiaations for the real geometry and for single-tone channel 
Bodel-wlth cross-sections of the averaged type. In the latter 
case - calculation differs from kinetic one so very little 
that it becoae possible to chart out a single coaaon line 
Instead of the two curves. Simultaneous ly, this common curve 
is vastly di/erent from the distribution not only in the lar- 
ger absorbtion region but in the moderator as well, prorided 
this distribution is representative of the fine channel struc- 
ture. In this way, the exact complex cell calculation should 
be based invariably on the detailed description of its geo- 
metry, this description on many occasions appearing more 
essential than the extensive description of diffusion pro- 
cess. fhe cases, however few they might be, when by force of 
accidental occurences the flux for a single zone model looks 
more or less coinciding with the actual channel flux (Pig. 5 ) 
the error in 0 might be substantial because of a suffi- 
cient neuron "depression" by the fuel layers. 

It is appropriate to mention that the approximation 
tiirns out to lack correctness so far as it is pwoorly descrij^ 
tive of the intra-channel neutron distribution, anyway, for 
the absorbing layers of the optical thickness. ( 2 ^) ^ 

Thus, the saw- like vibrations of the P^-flux in Pig . 3 
are about two times loss than the exact ones. Pig. 5 is par- 
ticularly dedicated (henceforth a little exaggerated repre- 
sentation) to the kinetic flux in the rings for the three va- 
lues ^ c of ring different from each other by two times; 
the flux alteration range in this case is almost strictly 
proportionate to Hq of ring. 

The comparison between the flux distributions calculated 
after the scheme under discussion wich the fluxes obtained 
in the suggestion of the neuti'on spectrum in che cell volume 

A ” n .7 ' r-, c CiTi f" r m o v'.w.i ^ i ^ ^ r C 1 \ 

, i_n ^ 

IS eiven lii Jigw 6 ^ ^ 




" 11 ^ 
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D«Bplt 0 the iacoaeidermbie divorgonce in fluxes the error 
In 0 looks eubotontially bxgi it lOf above all| due to the 
difference in the neutron gae temp-jraturee (See Table V). 

The fig* 7 is illuatrative of that the Pick law for the 
diffusion 3 provided the solution ia exact enough p 

does not hPId even qualitativeiy p and the nuaber of sign chan- 
ges of the neutron current ciurve J is not equal to the 
nunber of the extremuBS of the ^ flux* 

However, the diffusion current is by far closex* to the exact 
than Is the diffusion flux in relation to the exact flux, that 
can be helpful in the numerical solution for the calculation 
of derivatives in the formula (18). 

The flux redistribution in the process of the fuel bur- 
nout (Pig* 8) leads to a marked change in the absorption ra- 
tio in different rings. Th$s, if the heating power in the 
outer ring ouakes unity the middle and intter ring intensity 
changes with regard to the data of Table VI. 

Clearly, these redistributions must count in choosing 
the coolant flov rate through the gaps and the burnout calcu- 
lation simplification must be approached in a very careful 
manner. 

Pinally, Fig. 9-11 gives the neutron gas temperature 
distributions and those of neutron fluxes along; the cell ra- 
dius with two, three or four ring channels. One cannot fail 
to notice that the ”pr e-heating" depends 

very loosely on their number, and the spacing , for that mat- 
ter, and is primarily determined by the optical thickness of 
a ring. The calculated values of 0 for cell variants, as 
presented in the figures are brought together in Table VII • 

The inner fourth ring is so hardly screened that its ad- 
dition does not affect 9 , particularly at larger spacings 

when the thermal neutron generation in the Internal cavity 
does not matter. 

Table VII data were employed to analyse tlxe critical 
ezperiaents and proved rather benef itial for the agreement 
between the calculated and experimental data within the fra- 
mework of feir-group models. 
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1 

2 

3 

4 

5 

1.5 

0.91 40 

0.9228 

0.9229 

0.9245 

2.0 

0.8884 

0.8985 

0.8991 

0.9011 

3.0 

0.8389 

0.8512 

0.8529 

0.8554 

4.0 

0.7921 

0.8060 

0.8087 

0.8116 





TABLE V 


The Effect of Changing upon 0 

Pour Rings Two Rings 

"The 

Exact 


0.5741 


I Calculation” 



0.8499 

T 

n.g. 

= const. 


0.5969 

0.8570 





TABUS VI 


The Burnout and Changes 

of in 

the rings 

Average Uranium -235 

Middle 

Ring 

Inner Ring 

Burnout 




0 


0.5038 


O.3O8I 

25% 


0.5860 


0.3742 

50% 


0.6881 


0»4674 





TABLE VII 


Values 8 for Ring Channel Cells 


Cell 

Radius, Two Rings Three Rings 

Four Rings 


cm 





0.875 

0.57A 


0.879 

0 . 57 ^ 
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Comparison of the k-inetlc, 
r,- and F.-calcul ation of 
cell with ft three rin/^ 
channel (one-velocity 
npprox imation ) 


Comparison of the kinetic 
cell calculation results 
with a three rlnp: channel 
(Pif:,^) and the cell’model 
cftlculatlon 
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atSLizoxy units 
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O 1,0 to 9^0 kO O iJO kO kO 

Pig. 5 Pi^.6 


The effect of ring ‘’blackneae” The effect of neutron gas 
on the neutron flux(one-velo- temperature dlatributlon 
city approximation). Dotted on the equation solution 
line denotes the P--calcula- 
tlon of a slngle-zAne model 



The J current distribution The neutron flux changes 

in cells with two- and four- due to the three ring 
ring charmels(See the fluxes channel burnout 
in flg.6) 
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^ra«y wt/ij 
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